Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Zhang, T.*; 守田 幸路*; Liu, X.*; Liu, W.*; 神山 健司
Annals of Nuclear Energy, 179, p.109389_1 - 109389_10, 2022/12
被引用回数：1 パーセンタイル：29.26(Nuclear Science & Technology)The ID1 test was the final target test of the EAGLE experimental framework program. It was used to verify that during a core disruptive accident, the molten fuel could be discharged via wall failure of an inner duct in FAIDUS, a design concept for the sodium-cooled fast reactor. The ID1 results revealed that the wall failure behavior owed to the large heat flow from the surrounding fuel/steel mixture. The present study numerically investigated the heat transfer mechanisms in the test using the finite volume particle method in the three-dimensional domain. The thermal hydraulic behaviors during wall failure were reproduced reasonably. The present three-dimensional simulation mitigated inherent defects of our previous two-dimensional calculation and clarified that the solid fuel and liquid steel close to the outer surface of the duct can expose the duct to high thermal loads, resulting in the wall failure.
Zhang, T.*; 守田 幸路*; Liu, X.*; Liu, W.*; 神山 健司
Extended abstracts of the 2nd Asian Conference on Thermal Sciences (Internet), 2 Pages, 2021/10
日本のナトリウム冷却高速炉では、高速炉の炉心損傷事故における大規模炉心プール形成による再臨界を回避する方策として、内部ダクト付き燃料集合体(FAIDUS)が提案されている。本研究では、FAIDUSの有効性を実証するために実施されたEAGLE ID1炉内試験を対象に3次元粒子粒子法シミュレーションを行い、溶融燃料/スティールの混合プールからダクト壁への熱伝達機構を明らかにするための解析的検討を行った。
Zhang, T.*; 船越 寛司*; Liu, X.*; Liu, W.*; 守田 幸路*; 神山 健司
Annals of Nuclear Energy, 150, p.107856_1 - 107856_10, 2021/01
被引用回数：5 パーセンタイル：64.12(Nuclear Science & Technology)The EAGLE ID1 test was performed by the Japan Atomic Energy Agency to demonstrate the effectiveness of fuel discharge from a fuel subassembly with an inner duct structure. The experimental results suggested that the early duct wall failure observed in the test was initiated by high heat flux from the molten pool comprising liquid fuel and steel. In addition, the post-test analyses showed that the high heat flux may be enhanced effectively by molten steel in the pool. In this study, a series of thermal-hydraulic behaviors in the ID1 test was analyzed to investigate the mechanisms of molten pool-to-duct wall heat transfer using a fully Lagrangian approach based on the finite volume particle method. The present 2D particle-based simulation demonstrated that a large thermal load on the duct wall can be caused by direct contact of the liquid fuel with nuclear heat and high-temperature liquid steel.
吉田 啓之; 上澤 伸一郎
Proceedings of 2020 International Conference on Nuclear Engineering (ICONE 2020) (Internet), 7 Pages, 2020/08
The radioactive aerosol removal equipment is used as one of the safety systems of nuclear reactors. In this equipment, microparticles of aerosol are removed through gas-liquid interfaces of two-phase flow. The mechanism related to the removal of microparticles through the gas-liquid interface is not precise; a numerical evaluation method of performance of aerosol removal equipment is not realized. Then, we have started to construct a numerical simulation method to simulate the removal of microparticles through gas-liquid interfaces. In this simulation method, a detailed two-phase flow simulation code TPFIT is used as the basis of this method. TPFIT adopts an advanced interface tracking method and can simulate interface movement and deformation directly. Also, to simulate the movement of particles, the Lagrangian particle tracking method is incorporated. By combining the interface tracking method, and the Lagrangian particle tracking method, the interaction between interfaces and microparticles can be simulated in detail. To solve the Lagrangian equations of particles, fluid properties and fluid velocity surrounding aerosol particles are evaluated by considering the relative position of particles and gas-liquid interface, to simulate particle movement near the interface. In this paper, we show an outline and preliminary results of this simulation method.
橋本 慎太郎; 佐藤 達彦
Journal of Nuclear Science and Technology, 56(4), p.345 - 354, 2019/04
被引用回数：5 パーセンタイル：48.18(Nuclear Science & Technology)モンテカルロ法に基づいた粒子輸送シミュレーションは、加速器施設の遮蔽計算等の様々な目的で利用されている。モンテカルロ法による計算結果の信頼性を定量的に評価するためには、試行回数によって決まる統計的不確かさに加えて、計算で使われる反応断面積や計算の入力情報となる遮蔽材の密度がもつ誤差の影響を系統的不確かさとして求める必要がある。本研究の評価方法は分散分析に基づいており、これを鉄やコンクリートを遮蔽材とする中性子の遮蔽計算の解析に適用し、統計的及び系統的不確かさの両方を評価可能であることを示した。その際、入力情報の値を誤差の範囲で変動させるランダム条件法の他、三条件(中央値,上限値及び下限値)のみ変動させる三条件法を提案した。ランダム条件法は計算体系に関わらず適切に系統的不確かさが評価できるものの長い計算時間を必要とするのに対し、入力情報の誤差の影響が複雑な場合を除いて、三条件法は計算時間を抑えてランダム条件法と同じ評価結果を与えることがわかった。さらに、試行回数を増加させた場合の収束状況を判断できる新しい基準値を示し、必要最小限の計算時間で収束した評価結果が得られることを明らかにした。
守田 幸路*; 小川 竜聖*; 時岡 大海*; Liu, X.*; Liu, W.*; 神山 健司
Proceedings of 12th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics, Operation and Safety (NUTHOS-12) (USB Flash Drive), 11 Pages, 2018/10
EAGLE炉内ID1試験は日本原子力研究開発機構によって実施され、FAIDUSと称される内部ダクト付き燃料集合体からの早期燃料流出を模擬したものである。試験で生じた早期ダクト破損は、燃料とスティールから構成される溶融プールからの高い熱流束によるものと解釈されている。試験後の分析からは、壁面に燃料クラストが形成されない状況において、高い熱伝導度を有するプール中の溶融スティールによって溶融プールからダクトへの伝熱が効果的に促進されたことが示唆されている。本研究では、多成分多相流の粒子法に基づいた完全ラグランジェ法を用いて溶融プールからダクト壁への熱伝達機構を分析した。プール中の溶融スティールと燃料の混合と分離挙動およびこれらの挙動がプールからダクトへの伝熱に与える影響を調べるため、燃料ピンの崩壊、溶融プールの形成およびダクト壁の破損に至る一連の挙動を模擬した。現在の2次元粒子法シミュレーションでは、10MW/m
を超える壁面への大きな熱負荷は、核発熱を伴う液体燃料が壁面へ直接接触することによるものであることが示された。
Li, X.; 佐藤 一憲; 山路 哲史*; Duan, G.*
Proceedings of 26th International Conference on Nuclear Engineering (ICONE-26) (Internet), 8 Pages, 2018/07
溶融コリウム・コンクリート相互作用(MCCI)は軽水炉の仮想的シビアアクシデント時の後期フェーズにおいて炉容器外で生じる可能性のある重要事象である。本研究では、MPS法を用いてKITによって実施された模擬物質による成層化溶融プールの実験COMET-L3に対する3次元解析を行った。コリウム/クラスト/コンクリート間の伝熱は粒子間の熱伝導モデルで模擬した。さらに、ケイ酸系コンクリートではケイ酸系析出物の効果によって軸方向と径方向の浸食が異なる可能性が既往研究から示唆されていることから、2つの異なる解析ケースを実施した。解析の結果、MCCIにおいて金属コリウムは酸化物コリウムと全く異なるコンクリート浸食パターンを示しており、アクシデントマネジメントにおける格納系境界の溶融貫通時間の評価に考慮する必要があることが分かった。
Guo, L.*; 守田 幸路*; 飛田 吉春
Journal of Nuclear Science and Technology, 53(2), p.271 - 280, 2016/02
被引用回数：7 パーセンタイル：54.81(Nuclear Science & Technology)For the analysis of debris behavior in core disruptive accidents of liquid metal fast reactors, a hybrid computational tool was developed using the discrete element method (DEM) for calculation of solid particle dynamics and a multi-fluid model of a reactor safety analysis code, SIMMER-III, to reasonably simulate transient behavior of three-phase flows of gas-liquid-particle mixtures. A coupling numerical algorithm was developed to combine the DEM and fluid-dynamic calculations, which are based on an explicit and a semi-implicit method, respectively. The developed method was validated based on experiments of water-particle dam break and fluidized bed in systems of gas-liquid-particle flows. Reasonable agreements between the simulation results and experimental data demonstrate the validity of the present method for complicated three-phase flows with large amounts of solid particles.
n waves in tokamaksBierwage, A.; 篠原 孝司
Physics of Plasmas, 21(11), p.112116_1 - 112116_21, 2014/11
被引用回数：23 パーセンタイル：73.33(Physics, Fluids & Plasmas)The resonant interactions between shear Alfvn waves and resonant fast ions are analyzed numerically in realistic geometry and with a realistic particle distribution for a JT-60U plasma driven by two 400 keV N-NBs. In order to deal with the large magnetic drifts of the fast ions, a new orbit-based resonance analysis (ORA) method is developed and applied. The ORA method uses mappings that allows us to unify resonance conditions for trapped and passing particles, determine which harmonics are driven, and which orders of the resonance are involved. After analyzing the resonance conditions, the effects of wave-particle trapping, kinetic compression and linear resonance overlap are examined. Finally, we discuss the implications for frequency chirping, convective amplification, the long-time evolution of the system in the presence of a fast ion source, and for the interpretation of experimental observations.
滝塚 知典; 細川 哲成*
Contributions to Plasma Physics, 46(7-9), p.698 - 703, 2006/09
被引用回数：14 パーセンタイル：43.82(Physics, Fluids & Plasmas)ELMyHモードプラズマ中のELM崩壊後のダイバータ板への大きな粒子束と熱束は、トカマク核融合炉運転において重要問題である。このときのSOLとダイバータプラズマの過渡的挙動に及ぼす運動論的影響についてよくはわかっていない。ここで、先進的粒子シミュレーションコードPARASOLを用いてこの問題を研究する。粒子と熱の伝搬の衝突度への依存性を系統的に調べる。また、粒子リサイクリングの影響についても調べる。
松本 太郎; 内藤 裕志*; 徳田 伸二; 岸本 泰明*
Nuclear Fusion, 45(11), p.1264 - 1270, 2005/11
被引用回数：14 パーセンタイル：43.56(Physics, Fluids & Plasmas)核融合プラズマの中心部のような高温になると電気抵抗が小さくなるとともに、電子慣性等の非衝突効果(運動論的効果)が磁気再結合過程に深くかかわる。本研究では、これらプラズマの粒子性を物理的に忠実に取り扱う観点から、ジャイロ運動論的粒子モデルを用いた3次元シミュレーションにより、反転磁気シア配位トカマクにおける磁気再結合を伴うMHD不安定性に対する運動論的効果の解明を目的とした。反転磁気シア配位において二つの共鳴面の距離が近いと、互いの共鳴面における摂動が結合し、ダブルテアリングモード(DTM)と呼ばれる内部モードが生じる。高温プラズマにおいて電子慣性により誘起される運動論的DTMは、抵抗性モデルにより予期されるものより、速い内部崩壊を引き起こし得ることが明らかとなった。さらに、内部崩壊により掃き出された電流は、DTMが作り出した静電ポテンシャルに導かれて凹型に再集中し得ることが明らかになった。本研究は、低ベータ負磁気シア実験におけるディスラプション機構の一端を示しており、またDTMを経た電流分布の再分配により負磁気シア配位を再構築し得ることを示唆している。
寺田 宏明; 茅野 政道
Proceedings of 2nd International Conference on Radioactivity in the Environment, p.15 - 18, 2005/10
原子力緊急時において放射性物質の大気拡散のリアルタイム予測を行う計算モデルは非常に有用である。発表者らはこれまで緊急時環境線量情報予測システムSPEEDIとその世界版WSPEEDIを開発してきた。従来のWSPEEDIは質量保存風速場モデルWSYNOPと粒子拡散モデルGEARNから構成されていた。WSYNOPは診断型モデルであり気象場を予測することは不可能であるため、予測精度及び解像度が入力気象データに依存し、また鉛直拡散及び降雨沈着過程の詳細な考慮が不可能であった。これを改良するために大気力学モデルMM5を導入した。この改良版WSPEEDIを1986年に発生したチェルノブイリ原子力事故に適用し、数値モデルの予測性能の検証を行った。CASE-1:ヨーロッパ全域を含む広域計算と、CASE-2:広域とチェルノブイリ周辺域での2領域ネスティング計算の2ケースの計算を行った。CASE-1の計算結果より
Csの大気中濃度と地表沈着量を水平分布図及び統計解析で測定データと比較した結果、ヨーロッパスケールでの輸送の挙動をよく予測できていた。CASE-2の計算では、領域ネスティング計算によりCASE-1の広域計算では計算不可能であった詳細な沈着量分布を予測することができた。
松本 太郎; 内藤 裕志*; 徳田 伸二; 岸本 泰明
Physics of Plasmas, 12(9), p.092505_1 - 092505_7, 2005/09
被引用回数：2 パーセンタイル：6.85(Physics, Fluids & Plasmas)磁気流体力学(MHD)不安定性の研究は、プラズマを閉じ込める磁気面の保持/再生にかかわる、トカマク炉心プラズマ研究の中心課題である。核融合プラズマの中心部のような高温になると電気抵抗が小さくなるとともに、むしろ電子慣性,有限ラーマ半径等の非衝突効果(運動論的効果)が磁気再結合過程に深くかかわる。本研究では、ジャイロ運動論的粒子モデルを用いた3次元シミュレーションにより、磁気再結合を伴うMHD不安定性に対する運動論的効果の解明を目的とした。電子慣性の特徴的長さは電子スキン長であるが、それと比較して磁気再結合にかかわる他の素過程の特性長(イオンラーマ半径,抵抗層の厚さ等)が小さい場合、キンクモードは主として電子慣性によって誘起されると考えられる。しかし、このような電子慣性が支配的な運動論的キンクモードに対して、磁力線方向の電子圧縮性がわずかに存在する場合は、磁気再結合点近傍の電流構造を非線形的にY型からX型に変化させる。それゆえ、線形成長過程において電子慣性により支配されていたモードは、非線形的に成長率を加速させ、ある一定以上の成長率に達して内部崩壊に至ることが明らかとなった。また、このような非線形加速は、DTM等の内部モードに共通して現れる現象であることも判明した。
シミュレーション井戸村 泰宏
プラズマ・核融合学会誌, 81(8), p.581 - 592, 2005/08
ジャイロ運動論的粒子シミュレーションはトカマクプラズマの微視的乱流を研究するための有効な手段である。ジャイロ運動論的粒子シミュレーションにおける標準的な手法となっている法は、粒子ノイズの低減により粒子シミュレーションの効率を大幅に向上させ、完全トーラス配位の乱流シミュレーションを可能にした。本稿では法を概説し、完全トーラス配位のジャイロ運動論的粒子シミュレーションにおける問題点を議論する。
Cs deposition due to the Chernobyl nuclear accident寺田 宏明; 茅野 政道
Journal of Nuclear Science and Technology, 42(7), p.651 - 660, 2005/07
被引用回数：18 パーセンタイル：74.82(Nuclear Science & Technology)本研究の第1報(Terada 
2004)では、世界版緊急時環境線量情報予測システム「WSPEEDI」の改良版が大気力学モデルMM5とラグランジュ型粒子拡散モデルGEARN-newとのモデル結合により開発され、チェルノブイル原子力事故への適用計算によって大気中Cs濃度の予測性能について数値モデル(MM5/GEARN-new)の妥当性が示された。第2報である本論文では、このMM5/GEARN-newの予測性能評価をチェルノブイル原子力事故時のヨーロッパ域における降水量及び地表Cs沈着量の測定値を用いて行った。計算結果が測定値と比較された結果、MM5/GEARN-newは氷相過程を考慮した雲物理スキームにより高精度な降水量が予測された場合、ヨーロッパスケールでの沈着量分布を良い精度で予測することができた。さらに、ネスティング機能を用いてチェルノブイル周辺域についての高解像度計算が行われた。この結果、MM5/GEARN-newは従来の計算モデルでは計算できなかったチェルノブイル周辺での詳細なCs沈着量分布を現実的に予測することができた。
Villard, L.*; Angelino, P.*; Bottino, A.*; Allfrey, S. J.*; Hatzky, R.*; 井戸村 泰宏; Sauter, O.*; Tran, T. M.*
Plasma Physics and Controlled Fusion, 46(12B), p.B51 - B62, 2004/12
被引用回数：29 パーセンタイル：66.92(Physics, Fluids & Plasmas)ジャイロ運動論に基づく近年の第一原理プラズマ乱流シミュレーションの現状をレビューするとともに、ジャイロ運動論におけるv//非線形性,ジャイロ運動論的Poisson方程式における幾何効果の厳密な取り扱いといった、これまでの近似的に無視されてきた問題がシミュレーション結果に及ぼす影響を定量的に議論した。これらの効果の注意深い取り扱いの結果、エネルギー・密度保存を満足する乱流シミュレーションが可能になった。新しいシミュレーションでは(1)乱流の広がりが熱流速のバーストと圧力分布の局所平坦化の連鎖的な雪崩現象で生じる,(2)非線形の駆動されるE
B流はイオンラーマー半径の30倍程度のスケール長のグローバルなシア流になるといった、以前のシミュレーションでは見られなかった現象が明らかになった。
井戸村 泰宏; 徳田 伸二; 岸本 泰明
Journal of Plasma and Fusion Research SERIES, Vol.6, p.17 - 72, 2004/00
トカマクプラズマにおいてITG-TEM(イオン温度勾配駆動-捕捉電子モード)乱流により生じるイオン及び電子異常輸送を調べるためにグローバルジャイロ運動論的トロイダル粒子コードGT3Dを開発した。このコードによる予備的なITG-TEMの線形計算においてITG-TEMの基本的な性質を確認した。捕捉電子はITGモードの成長率を増大させるだけでなく、別の電子モード、TEMを励起する。
あるいは
により支配的なモードはITGからTEMに変化する。サイクロンパラメータを用いた線形ベンチマークにおいてGT3D, GTC(PPPL-CUI)、及びFULL(PPPL)コードの間で良い一致が得られた。
O
Mo, S.-K.*; Denlinger, J. D.*; Kim, H.-D.*; Park, J.-H.*; Allen, J. W.*; 関山 明*; 山崎 篤志*; 角野 宏治*; 菅 滋正*; 斎藤 祐児; et al.
Physical Review Letters, 90(18), p.186403_1 - 186403_4, 2003/05
被引用回数：143 パーセンタイル：95.48(Physics, Multidisciplinary)複雑な相図を示すVOについて、その金属相での光電子スペクトルにおいて顕著な準粒子ピークをはじめて観測した。また、この明瞭なピークが温度効果とともに成長することを示すために、局所密度近似と量子モンテカルロシミュレーションを用いた動的平均場理論を組み合わせた新しいスペクトル計算結果を報告する。実験で得られたピーク幅,ウェイトはともに理論より相当大きい。
井戸村 泰宏; 徳田 伸二; 岸本 泰明
Nuclear Fusion, 43(4), p.234 - 243, 2003/04
被引用回数：119 パーセンタイル：94.96(Physics, Fluids & Plasmas)ジャイロ運動論的トロイダル粒子コードGT3Dを開発した。GT3Dでは、軸対称トロイダル系における運動の恒量により定義される正準マックスウェル分布に基づく手法,高
モードのグローバル解析を可能にする準バルーニング展開、あるいは、粒子,エネルギーなどの保存則を改善する最適化粒子配分法といった新しい手法が実装されており、実装置パラメータにおける高精度のプラズマ乱流シミュレーションが可能になっている。本コードを大型トカマク
におけるITG(イオン温度勾配駆動)モードの解析に適用し、反転磁気シア配位は
領域におけるITGモードの実質的な安定化に寄与することを明らかにした。また、GT3Dでは正準マックスウェル分布を用いた実装により、軸対称モード、特に、乱流抑制に寄与する乱流駆動
帯状流の応答が正しく表現されている事を確認し、従来の局所マックスウェル分布を用いた乱流シミュレーションにおける重要な問題点を明らかにした。
木田 孝; 梅田 幹; 杉川 進
JAERI-Data/Code 2003-001, 29 Pages, 2003/03
JAERI-Data-Code-2003-001.pdf:6.34MB
燃料サイクル安全工学研究施設(NUCEF)では、臨界実験で使用する硝酸プルトニウム溶液を調製するために、銀媒体電解酸化法によるウラン・プルトニウム混合酸化物(MOX)粉末の溶解を計画している。この溶解運転の支援のために、MOX溶解シミュレーションコードの開発を行なった。本報は、開発したシミュレーションコードの概要、実験データとの比較及び溶解速度に及ぼすMOX粉末溶解パラメータについてまとめたものである。本コードは、Ag
によるPuO粉末溶解に対するZundelevichのモデルに基づいて作成し、Agの収支式に亜硝酸による影響を加えるとともに、MOX粉末の表面積は粒径分布により算出する方式を採用した。また、コードの妥当性を確認するために、実験データとの比較を行なった。この結果、本コードにおいて溶解速度定数等のパラメータに適切な値を用いることで、溶解挙動を良く模擬できることを確認した。また、溶解速度に及ぼす影響が最も大きいパラメータは、MOX粉末粒径であることが判明した。
